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Abstract 



^D . Studies on nanoscale materials merit careful development of an electrostatics model concerning 

H . discrete point charges within dielectrics. The discrete charge dielectric model treats three unique 

interaction types derived from an external source: Coulomb repulsion among point charges, direct 
polarization between point charges and their associated surface charge elements, and indirect po- 
larization between point charges and surface charge elements formed by other point charges. The 
model yields the potential energy, U{N), stored in a general A^ point charge system differing from 



C^ • conventional integral formulations, 1/2 J E ■ DdV and 1/2 j p^dV , in a manner significant to the 



treatment of few electron systems. 



I. INTRODUCTION 

The most common conceptual approach to evaluate the total electrostatic potential en- 
ergy stored in a dielectric containing free point charges is the method of charge assembly. 
Textbook treatments^"- replace discrete point charges within a dielectric with a continuous 
free charge density, p. Energy stored in a dielectric containing p is easily evaluated by in- 
tegration over the volume of the dielectric in which the smallest increment of charge, dq, 
is zero. Though this is a reasonable approximation when p consists of a large number of 
point charges," it fails to reflect the physical nature of point charge configurations commonly 
found in nanoscale systems today. Moreover, the imposition of a continuous charge density 
prohibits the practical formation of any realistic charge distributions composed of infinitesi- 
mal, non-existent charge increments. The present model is set apart from these approaches 
in that a summation of interaction terms is obtained rather than integral approximations 
based on the interaction of p with a macroscopic total potential, $,- 

U{p) = \j^pMV (1) 

or an electric field formulation,- 

UiE) = - f E- DdV. (2) 

2 Jv 

The objective of the present model is to account for the interaction of each new point 
charge to the system with the sum of all electrostatic potentials due to every pre-existing 
charge and charges formed during the assembly process without assumption of an infinitesi- 
mal charge increment, dq, to refiect the natural existence of a fundamental charge increment, 
e. 

The discrete charge dielectric model has been used for numerical evaluations in re- 
cent studies concerning the capacitance of few-electron spherical quantum dots^^i""— and 
"meta-atoms" ^-, cited as leading to charge carrier localization near quantum dot surfaces^ 
and played a central role in the discovery of a classical electrostatic fingerprint of atomic 



structure.— li^ The present model represents a redevelopment and generalization of a previ- 
ously published interactions picture initially used to evaluate the electrostatic energy of 
a silicon sphere containing one or two electrons embedded in a uniform silicon dioxide 
environment.— li^ Exclusion of one-half of the "self-polarization" interaction energy term 
in the original picture is shown to be inconsistent with the current formulation in which 
a more appropriate description, "direct polarization," is used to both eliminate confusion 
with self-interaction or self energy - interaction of a point charge with its own potential - 
and to differentiate it from indirect polarization interactions ( "polarizations" in the original 
picture) which are not directly involved in the polarization of the dielectric. 

A practical definition of stored energy is that delivered to the system by an external 
agent acting on point charges during assembly and must be equal to that recovered by an 
external conservative energy source when the system is dismantled. This is necessary such 
that the dielectric is returned to its original, neutral state. The amount of energy displaced 
by polarization of atomic or molecular constituents of the dielectric must likewise be restored 
to its internal conservative energy source (atoms or molecules) but does not contribute to 
the energy stored in the system as it already exists within it. 

In the following, the electrostatic interaction between two point charges of like sign is 
carefully described. Seemingly trivial, the development is exceedingly instructive toward a 
physical intuition and constructs a framework in which the other types of interaction are 
developed. This approach identifies the mutual distribution of interaction energy between 
charges and categorizes interactions as directly resulting from work done by an external 
agent or a response of charges already in the system to new charge arrivals. 

The electrostatic constraint of a single point charge to a dielectric object is examined. 
Formation of a bound surface charge element at the dielectric interface is the result of 
direct polarization, an example of the charged system's response to a newly introduced 
point charge. In the process of introducing a second point charge to the system, the point 
charges interact directly with each other, and a new bound surface charge is formed by direct 
polarization. Additionally, the external agent must do work on the new point charge against 
the electrostatic potential due to the bound surface charge element formed exclusively by 
the first point charge. This represents an indirect polarization in which energy is mutually 



shared between the two point charges as the first point charge is hkewise subjected to the 
potential due to the additional bound surface charge element formed by the second point 
charge. The latter interaction is a necessary response to the formation of an additional 
bound charge. The energy associated with this response can only derive from work done 
by the external agent. For completeness, interactions among bound charges are arguably 
intrinsic to the system and have no direct source of energy as may only be supplied by the 
external agent. These interactions are therefore explicitly taken into account within the 
framework of the discrete charge dielectric model during assembly. 

For few-electron systems commonly found in nanoscale studies the generalized summation 
expression differs from Eqs. [T] and [2] by up to a factor of 1/2 (A^ = 1). This result is 
not unexpected as it is consistent with a similar charge-by- charge derivations^ of the so- 
called "quantum" capacitance expression, Cq^ant = Q"^ /U ■, which differs from the familiar 
"classical" capacitance expression, C^jg^gg = Q^/2f/, by a factor of 1/2. Similarly, a 1/2 
capacitance disparity between "open" (isolated) and "closed" (connected) single electron 
devices is known.— The sources of these disparities are discussed and represent one-half the 
total direct polarization energy in the present model. 

In the other extreme, as A^ grows large, the direct polarization energy term becomes a 
negligible fraction of the total energy as required for the approximation of dq — )■ 0. Hence, the 
discrete charge dielectric model is applicable to nanoscale systems in which the behavior of 
few free electrons dominates the electronic properties of a material, as well as many-electron 
systems limited by computing capacity and other practical considerations. The generalized 
expression may be numerically minimized to obtain a final electrostatic configuration for 
any A^-electron system in satisfaction of necessary boundary conditions and the shape and 
size of the dielectric. 

II. DISCRETE CHARGE DIELECTRIC MODEL 

The electrostatic force acting on a point charge, g, by an electric field, E, due to all other 
sources of charge in a system, is given by F = gE. The work done to displace the point 
charge by a small distance, dx, is dW = — F ■ dx such that the total work done to move q 



from a non- interacting location, r„, to an interacting location, r^, in the vicinity of a charged 
system is given by 

/"T" ■ PIT ' 

W = I dW = -q IE-d:s. 

Jrn Jr„ 

= q [$(r,) - <l>(r„)] 

= g$(n). (3) 

Here, $ is the electrostatic potential due to all sources of charge present in the system before 
and during the displacement of q. 

A single point charge in the absence of other charges is considered to have no stored 
energy as no work is done to move the point charge around. A point charge, qi, interacting 
with its own electrostatic potential. 
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at its location, Tj, is infinite in magnitude. This self-interaction energy is excluded from 
evaluation of the total stored energy as the external agent is not involved. 

A. Coulomb Interaction 

To evaluate the "Coulomb interaction"— energy between gi and q2 in free space having 
the same sign consider an external agent acting on q2 in the direction of gi . As g2 approaches 
gi retreats to infinity as there are no physical constraints acting on it. If there are physical 
constraints, an evaluation of the energy associated with the constraining forces must be 
made. This concept is central to the discrete charge dielectric model, and dismisses the 
unnecessary prescription of an initially "fixed" charge configuration conventionally implied 
in the electrostatics literature. Hence, the response of the pre-existing charge configuration 
to new charge arrivals is paramount. This perspective readily facilitates the introduction of 
new charges to the system. 



If gi is found within a dielectric characterized by dielectric constant, e, the displacement 
of gi is the result of gi interacting with the potential due to q2- This interaction represents 
the response of gi to the arrival of g2, but may only derive its energy indirectly from the 
work done on g2 by the external agent. As the interaction terms are equal by symmetry the 
energy supplied by the external agent acting on g2 against 0(gi,r2), 



Uciqi,Q2) = g20(gi,r2) 

= *T^^^. (4) 

ATTEoe In - r2| 

is mutually shared between gi and g2. The notation (/)(gi, r2) denotes the source of electro- 
static potential and location at which it acts on g2. The notation Uc{qi,q2) denotes the 
charges with which energy is mutually shared. 

The Coulomb interaction is subdivisible into two interactions: 1) a response of gi indi- 
rectly supplied energy through 2) direct Coulomb interaction of an external agent acting 
upon g2 during the assembly process. The latter is the sole source of stored energy involved 
in the Coulomb interaction. 

Fig. [1^ schematically depicts a configuration of two point charges and resulting dipole 
distributions of net-polarized bound surface charge elements, cxi and (J2, in a dielectric. 
Fig. [T)d shows the electrostatic interactions necessary for a complete evaluation of the total 
stored energy. The Coulomb interaction, Uc, is schematically shown with mutually-directed 
arrows between point charges gi and g2. 

B. Direct Polarization 

Consider the process by which gi is introduced into the dielectric in the previous discus- 
sion. The dielectric is initially electrically neutral as the positive and negative charges within 
the material are arranged such that there is no net accummulation of charge. To this system 
is introduced a negative point charge, gi, by an external agent. The external agent must 
supply energy in its direct interaction with gi to act indirectly against the restorative nature 
of the charged atomic or molecular constituents of the dielectric in a process of polarization. 



Due to the dielectric disparity at the boundary, a net negative region of charge is formed. 
The result is the formation of a unique bound surface charge element, ai. If gi were removed 
from the system, the restorative nature of the atomic or molecular constituents act to return 
the displaced energy to its internal source. This process is representative of a conservative 
internal energy source. 

The direct polarization energy, 

^dir(^i' ^i) = ?i0(^i) '^i)^ (5) 

is mutually shared between gi and cti, contributing to the energy stored in the free point 
charge distribution, Uf^.^^, and the bound surface charge, f^bound' respectively. This is 
schematically shown in Fig. [T|d with mutually-directed arrows for ^/^ir- "^^^ total energy 
stored in a single point charge system is given solely by the direct polarization interaction, 

t/(l) = t/dir(gi,^i). (6) 

C. Indirect Polarization 

Introduction of a second point charge, q2, to the system (Fig. [T^) involves work done on 
q2 by an external agent against the electrostatic potential of gi, o"2 and o"i. Interactions with 
potential due to gi and (J2 are respectively given by Eq. 1^ and 

f^dir(92, (^2) = g20(c^2, ^2). (7) 

Interaction with potential due to o"i results in the indirect polarization of o"i by g2 given by 

^ind(5i' ^2) = q2'P{(Ti,r2). (8) 

The bound surface charge element ai is exclusively dependent upon the magnitude and 
location of gi. Therefore, interaction of g2 with 0(cri,r2) does not directly redistribute ai. 



Instead, the introduction of g2 results in the formation of a^ which presents a potential that 
displaces q\. The displacement of q\ results in a redistribution of a\ by direct polarization. 
Hence, a\ is indirectly polarized by g2- 

Energy associated with the interaction of q\ with 0(o"2,ri) must have a source. Work 
done on g2 is the sole source of energy. The source of energy associated with the indirect 
polarization of 02 by q\ is found by distributing energy among the entire system of charges. 
Tracing through Fig. [TJd, the interaction energy associated with g2 acting against the poten- 
tial due to o\ remains associated with g2 (hence, the arrow from a\ to q-i). This energy may 
be distributed along two interaction "paths": 1) through Coulomb interaction with gi, and 
2) through direct polarization with o^- Energy associated with the indirect polarization of 02 
by gi remains associated with gi, explicitly denoted in Eq. [HI as a^ is exclusively dependent 
on g2- 

In satisfaction of interaction symmetry, energy is redistributed throughout the charge 
system (Fig. [Do) such that energy resulting from the indirect polarization of o\ by g2 is 
mutually shared by q\ and g2- The total indirect polarization energy contributes to t^free- 
Any contributions made to f^bound ^^^ taken into full account through respective direct 
polarizations due to mutual displacements of q\ and g2 by a^ and cri, respectively. 

The energy stored in a dielectric containing two point charges is then 

2 

U{2) = Uciqu q2) + U^ndi^i, 92) + E UdiriQh ^i)- (9) 

D. Surface Charge Elements 

For completeness, electrostatic interactions among surface charge elements, ai, have no 
direct source of energy and must instead draw upon energy provided indirectly by the respec- 
tive direct polarizations as-evaluated. Noted previously, the magnitude and distribution of 
each (Tj is exclusively dependent upon the magnitude and location of its associated g^. Simi- 
larly, energy associated with the indirect polarization interaction gi0(o"2, ri) is not evaluated 
directly but draws mutually upon the indirect polarization interaction g20(o"i,r2). Addi- 
tionally, careful treatment of Coulomb interaction energy between two point charges yields 

8 



only one-half of the apparent sum of energy resulting from both interactions. Explicitly 

Uc{qi, 52) = - [gi0(g2, n) + g20(gi, ^2)] (10) 

is equivalent to Eq. HI Eq. H] implicitly takes into account the interaction gi0(g2,?"i) as 
(i2<P{<li-,^2) is an indirect internal response. 

Hence, interactions among ai are explicitly taken into account by evaluation of all work 
done by an external agent acting upon charges assembled to the dielectric within the frame- 
work of the discrete charge dielectric model. 

E. N Electrons 

As each new electron is assembled to the system its interaction with the electrostatic 
potential due to all charges are taken into full account by three types of interaction. The 
model may be generalized to a dielectric containing A^ electrons. The total stored electro- 
static energy is given by generalization of Eq. |9l 



Y N N 1 Af TV N 

U{N) = ,:T.i: UciQ^, Qj) + :^ E E f^md(9^, 9.) + E f^dir(?M ^^) (H) 

where mutual energy distribution with charge pairs within each interaction is explicitly 
denoted. 



III. DISCUSSION 

The total energy stored in a dielectric containing discrete point charges must be equal to 
the amount of energy delivered to the system in the form of work done by an external agent 
during charge assembly. This energy is evaluated by taking into account the interaction 
of each discrete point charge with the total electrostatic potential arising from all charges 
in the system both before and during the assembly process. Full partitioning of energy 
during assembly is provided by dividing interactions into three types: Coulomb interactions 



among point charges, direct polarization between each point charge and its exclusively bound 
surface charge element, and indirect polarization resulting from the interaction of each new 
point charge with bound surface charge elements formed by other point charges. 

Subdividing each type of interaction into those contributing to the energy stored due 
to work directly done by an external agent and those that must be indirectly satisfied by 
response to new charge arrivals, partitioning of energy is completely comprehended. By 
the symmetric relationship between interacting charges, mutual distribution of energy is 
observed. The energy may be broadly associated with the free point charge distribution, 
[/fj.gg, and the bound surface charge, t^bound-" Energy resulting from Coulomb interactions, 
indirect polarization interactions and one-half of direct polarization interactions contribute 
to f/fj-ee- O'^^ h^lf of t^^ energy resulting from direct polarization interactions contribute 
to f^bound- '^^^ total bound charge, a, is comprised of contributions from bound surface 
charge elements, cxj. Interactions among ai are explicitly taken into account in the present 
model, having acknowledged them as internal responses with no direct sources of energy 
beyond those discussed. 

Unique to the discrete charge dielectric model is its treatment of discrete point charges 
in dielectrics leading to a general summation, Eq. [TTl in place of conventional integral 
formulations, Eqs. [1] and [21 The difference between the present formulation and integral 
formulations is the factor of 1/2 in the latter. The last term of Eq. [TT] is not halved as 
both the energy required to form surface charge elements and the energy associated with 
point charges due to interaction with resulting surface charge elements must be added to 
the system. 

The amount of work done in the formation of each a^ is given by^ 

W= (f){a,,n)dq= (12) 

Jo 2 

since (j){ai,ri) oc q, and ai is justifiably formed by charge increments dq — ?■ 0. One-half 
of each direct polarization interaction qi(j){ai, Tj) is consumed in the process and stored in 
ai, while the remaining half is associated with each g, (part of t^fj-ee)- The bound energy, 
f/bound' ^^ust not be removed by the invention of a negative sign on a "spring energy" term, 
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^sprin^ nor merely set to zero except in the restrictive case of fixed external source charges 
of an electric field.- For this reason there is no factor of 1/2 on the third term of Eq. [TTl 

In the few-electron extreme of A^ = 1, the discrete charge dielectric model yields twice 
the total stored energy given by textbook formulations. This result is anticipated as it leads 
to an expression of capacitance^"— identical to that proposed in the past couple decades 
for a so-called "quantum" capacitance expression, Cqua^t = Q"^ /U that differs from the 
conventional expression by a factor of 1/2. An explanation of this disparity is immediately 
identified with the presumption of (ig — )■ in the derivation of C'j.jg^gg given by evaluation of 
the work done to "charge" a dielectric, 



W =JdW=f (j)dq 
Jo 

2C 



/o^ ¥9= ^- (13) 



Alternatively, derivation using a discrete charge element, Ag = e leads to W = Q'^/C 
with Q = Ne using a charge-by- charge derivation similar to that leading to a "quantum" 
capacitance^^. 

Additionally, the one-half disparity between the present model and conventional formal- 
ism is also observed between isolated and connected single electron devices.— In this context, 
an imposed condition of fixed external source charges of the electric field in Eq. |2f is justified 
within a "closed" (connected) device as metallic electrodes rigidly define equipotential sur- 
faces from which capacitance is derived. Moreover, electrodes consist of many free electrons, 
though with today's nanowires and the like, presumption of fixed external source charges is 
no longer merited. Absent metallic electrodes, as with biological entities for instance, for an 
"open" (isolated) device, images of internal electrons play the role of external source charges. 
Consequently, energy associated with image displacement must be taken into account - the 
equivalent of one-half the total direct polarization energy. This constitutes fully one-half 
the energy stored in an isolated single electron device. 

As A^ — )■ oo, the direct polarization term approaches a negligible fraction of the total 
stored energy as justified by the approximation dq — )■ 0. The present model applies to a 
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general A^ electron system. Its merit is significant to the treatment of electrostatic properties 
of few-electron devices. 

The summation formulation provided by the discrete charge dielectric model is amenable 
to numerical solution. The total stored energy of a system, Eq. [TTl is minimized to obtain 
electrons in their final electrostatic locations. Used in the treatment of spherical quantum 
dots, the discrete charge dielectric model has played a centrol role in the identification of a 
classical electrostatic "fingerprint" of the periodic table of elements.— li^ The "fingerprint" 
consists of a unique non-uniform distribution of electrostatic potential energy intrinsically 
associated with point-charge symmetries of neighboring iV-electron systems subject to spher- 
ical constraints (e.g. a central field or a spherical quantum dot). These non- uniformities are 
consistent with changes in inter-atomic subshell orbital geometries, such as between the 
spherical s-orbital of an outer electron in beryllium {Z = 4) and the energy associated with 
the dumbbell-shaped p-orbital of neighboring boron {Z = 5) as may be empirically char- 
acterized by their respective scale-independent ionization energies. Outside of this work, 
classical electrostatic shell-filling behavior has not been observed in treatments of spherical 
quantum dots, or "artificial atoms," though at least one unsuccessful attempt is described 
in the literature.— 
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FIG. 1. Discrete Charge Dielectric Model a) Two point charges, qi and q2, repel each other and 
polarize the dielectric such that net bound surface charge elements, 0"i and (72, are formed, b) 
Coulomb repulsion, Uc, is schematically represented by a mutually-directed arrow between qi and 
(72 indicating mutual energy distribution. Likewise, for direct polarizations, C^fjjp, energy is mutually 
shared between each point charge, qi, and its respectively polarized surface charge element, qi. As 
(7j are exlusively dependent upon respective qi energy derived from indirect polarizations, U^j^^, 
are mutually shared between qi and q2, denoted by singly-directed arrows. 
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